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Abstract 

We present a new method for multi-component power spectra estimation in multi- 
frequency observations of the CMB. Our method is based on matching the cross and 
auto power spectra of observation maps to their expected values. All the component 
power spectra are estimated, as well as their mixing matrix. Noise power spectra 
are also estimated. The method has been applied to full-sky Planck simulations 
containing five components and white noise. The beam smoothing effect is taken 
into account. 
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1 Introduction 

Precise measurements of Cosmic Microwave Background (CMB) angular power 
spectrum, is one of the main objectives of modern Cosmology. Several experi- 
ments, such as ACBAR [10], Archeops[2], Boomerang [11], CBI [13], DASI [14], 
MAXIMA [12], VSA [15] and more recently WMAP [16] already give precise 
results over a large range of angular scales. An unbiased estimate of the CMB 
angular power spectrum requires important efforts in data processing because 
of the presence of several sources of contamination in the observations. Fore- 
grounds, in particular (Galactic dust, synchrotron, extragalactic point sources 
and Sunyaev-Zel'dovich effects), will be the main contaminants for CMB mea- 
surements, such as that of the future Planck mission. Fortunately, the avaibii- 
ity of several detectors operating at several frequencies, ranging from 30 GHz 
to 857 GHz for Planck, permit the separation of most astrophysical compo- 
nents due to their different emission laws. 

Starting from the different observation maps, current methods for CMB power 
spectrum estimation consist first in making the cleanest CMB anisotropy map 
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possible, using component separation techniques, like Wiener filtering [1,9,7] 
or MEM [5,8]. Spectral estimation is then performed on the individual com- 
ponent maps. This approach is not fully satisfactory for two reasons: First, 
component separation requires prior knowledge of the electromagnetic spectra 
of the components, which are not necessarily very well known. Second, after 
component separation, in the subsequent spectral estimation step, it is nec- 
essary to remove the power spectrum of the residual noise in the extracted 
CMB map. A small error in the noise evaluation introduces a bias in CMB 
power spectrum estimation. 

A new approach has been considered in paper [4]. This approach is based on 
the analysis of cross- and auto-power spectra of different observation maps by 
different detectors. Considering the simple case where only CMB anisotropies 
are present with noise, the cross-power spectra give a direct measurement 
of the CMB power spectrum, assuming that the noise is uncorrelated be- 
tween detectors. In this limited context, the technique has been used by the 
WMAP team. Our approach deals simultanously with multi-component data. 
By jointly analysing the different observation maps, we estimate the power 
spectra of all the components, without making any priors on their emission 
laws. Thus, the contribution of each component to each detector, given by 
the mixing matrix A (e.g. (1)), is estimated directly from the data. The noise 
power spectra are also estimated. The method is based on likelihood maxi- 
mization in the Whittle approximation assuming that the observations are a 
linear mixing of independant components and independant noise. It can be 
viewed as a blind multi-detector multi-component (MDMC) spectral match- 
ing, in which the spectral diversity of the various components is used. The 
method has been implemented on harmonic coefficients of all sky maps. We 
account for the finite spatial resolution of the detectors. In this paper, we 
describe the basics of the spectral matching method and we present its appli- 
cation on full-sky simulated Planck observations. 



2 Model for the observed sky 



The key assumptions are that the sky emission at a given frequency is a linear 
superposition of astrophysical components, and that their emission laws do 
not depend on sky position. The signal measured by a detector is the sky 
emission convolved by a beam shape (depending on the detector), plus an 
additive noise. A spherical harmonic expansion over a full sky map, observed 
by detector d, gives us: 

x d (£, m) = b d (£) A dc s c (£, m) + n d (£, m) (1) 

c=l 
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where s c is the emission template for source c, rid represents the noise and A is 
the mixing matrix. Each element of the mixing matrix results from the integra- 
tion of the emission law of one component over one detector frequency band. 
The coefficient bd{£) refers to a Legendre Polynomial expansion of the beam, 
depending only on r.r' (we assume a symmetric beam). Let us consider the di- 
agonal matrix B(£) such that the diagonal element Bdd{£) = bd{£). It is useful 
to define the "deconvolved" observation coefficients x'(£,m) = B(£)~ 1 x(£,m), 
that we write in matrix form: 

x'(£,m) = As(£,m) + B(£)- 1 n(£,m) (2) 



Observation power spectra 

Our method is based on the analysis of cross- and auto-power spectra of the 
"deconvolved" observation maps, whose ensemble averages are expressed as : 

R x ,{£) = £ (x'(£,m)x'(£,my). (3) 

' m=—t 



where denotes transpose-conjugation. We define C{£) and N{£) as the com- 
ponent and the noise power spectra respectively, we assume statistical inde- 
pendence between components and also between the noise of different detec- 
tors, so C{£) and N(£) are diagonal matrices. The expected observation power 
spectra become, using (2): 

R x ,(£) = AC(£)A l + M(£) (4) 



where M(£) = B{£) 2 N(£) is also a diagonal matrix. 
Observation power spectra are estimated in the data by: 



= £ x'(£,m)x'(£,my (5) 

m=—i 



Expected and measured power spectra can be averaged over bins q such that, 
for example, £ m i n (g) < £ < £ max (q). Their structures, as in equations (4) and 
(5) are preserved: 



R x ,(q)=AC(q)A t + M(q) (6) 

1 ^max(lj) t 

Rx'(q) = — E E x'(£,m)x'(£,my (7) 
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where n q is the number of modes in bin q. According to equation (6), the 
non-diagonal part of R x >(q) contains only contributions from components (by 
supposition of noise independance) . This shows that the cross power spec- 
tra gives a direct estimate of the mixing of component power spectra; this 
mixing can be inverted, even without priors on A, assuming independance of 
components and some spectral diversity. 



3 Multi-component spectral matching 

Our method is based on minimizing the mismatch between the empirical power 
spectra of the observations (7) and their expected values (6), in order to mea- 
sure simultenaously the averaged component power spectra C(q) and the mix- 
ing matrix A (see discussion in paper [6] on a blind vs a semi-blind approach) . 
As we wish to avoid using any priors, the averaged noise power spectra M(q) 
are also estimated. These parameters are referred to as 9 = {A, C(q), M(q)}. 
The mismatch is quantified by the average divergence measure between the 
two matrices: 

<$>{9) = Y.n q D{RAq)iRA<l)) (8) 

q=l 

Assuming that the spherical harmonic coefficients of the components are ran- 
dom realizations of a Gaussian field of variance R x /(q), and that they are all 
uncorrelated, the log-likelihood (up to an irrelevant factor) takes the same 
form as in equation (8) (in the frame of the Whittle approximation), in this 
case the divergence is given by D(R 1 , R 2 ) = tr(i? 1 i? 2 _1 ) — log det(i?i-R^ 1 ) — m d 
[4]. The estimate 9 is such that <!>(#) is minimum. The connection with the 
log-likelihood garanties asymptotically minimum variance of the estimates and 
the absence of bias (this last property is preserved even if the components and 
noise are non-stationnary and non-Gaussian). 

In practice, we optimize the criteria <£>(#) in a first step using an EM (Expectation- 
Maximization) algorithm and the convergence is ended in a second step using 
a Quasi- Newton method (see papers [4,3] for a description of the algorithm). 



4 Application to Planck observation simulations 

We now turn to the application of MDMC spectral matching method on sim- 
ulated all-sky Planck observations provided by the Planck consortium. The 
maps are generated for the ten frequency channels of the Planck instruments. 
Five components : CMB, thermal dust, the two SZ effects and synchrotron, 
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and white noise at the nominal instrument level are mixed according to equa- 
tion (1). The beam sizes are also taken at their nominal values. See paper [8] 
for more details on the simulations. 

The cross- and auto-power spectra of simulated Planck observations are com- 
puted following equation (7), up to the multipole t = 3000 and using bins of 
width Al = 10. We try to estimate four components. Indeed, the kinetic SZ 
can not be separated from the CMB anisotropies with our method since these 
two components have proportional emission laws (they form one component). 
We estimate all component power spectra and all mixing matrix elements ex- 
cept three, corresponding to three of the four elements at 857GHz, which are 
fixed at zero (we expect only thermal dust at this frequency). This operation 
does not affect CMB power spectrum estimation results, but allows us to break 
degeneracies between galactic components having almost proportional power 
spectra. We assume as prior information that the noise is white (N indepen- 
dant of I). The total number of parameters is 10 x 4 - 3 + 4 x 300 + 10 = 1247, 
compared to 300 x 10 x (10 + l)/2 = 16500 data elements. 
The mixing parameters relative to the four components are estimated with 
good accuracy (see paper [6]). In particular, we are able to put strong con- 
strains on emission laws of Galactic components. Figure 1 shows the es- 
timated CMB power spectrum and the estimated relative errors given by 
(\C(q) — C(q)\/C(q)). It appears that the method allows us to accurately 
estimate the CMB power spectrum up to £ ~ 2000. No bias in the estimation 
is seen over the entire I multipole range. 



5 Discussion 
Single component 

The sensitivity of current CMB experiments lead us to expect that at very 
high galactic latitude at frequencies around 150 GHz, other astrophysical com- 
ponents are negligible as compared to CMB and noise (this will not be the case 
for Planck). Even in this simple case, our method, applied on small sky re- 
gions observed by multi-detector experiments, gives comparable or even better 
CMB power spectrum estimates than the "classical" techniques. The advan- 
tage of our method is that an independant noise power spectra estimation is 
not required since we explicitly consider cross power spectra of observation 
maps. 

Component map separation 

Our spectral matching method yields all the parameters needed to implement 
a Wiener-based or MEM component separation on the maps. Note that we 
proceed in making a spectral estimation followed by a component separation, 
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Fig. 1. Left: Estimated CMB power spectrum and CMB power spectrum of the 
input map. Right: Relative errors made on the CMB power spectrum estimation. 

instead of the inverse, as curently done in classical approaches. We have ap- 
plied Wiener filtering on simulated Planck observations using the estimated 
parameters. Component maps are accurately estimated, in particular output 
CMB map does not show residual contaminations in galactic plane regions. 



6 Conclusion 

We have presented our blind spectral matching method for all-sky multi- 
detector CMB observations. The main objective is to measure the power spec- 
tra of all the components, including the CMB, and their contribution levels 
at each observation frequency. The method exploits in particular the cross- 
power spectra of observation maps, giving an unbiased estimate of component 
power spectra. The method has been applied on full-sky Planck simulations 
containing five components and white noise. The power spectrum of the CMB 
is accurately estimated up to I ~ 2000 in bins of size Al = 10. 
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